Lattice distortions in a sawtooth chain with Heisenberg and Ising bonds. 
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An exactly solvable model of the sawtooth chain with Ising and Heisenberg bonds and with cou- 
pling to lattice distortion for Heisenberg bonds is considered in the magnetic field. Using the direct 
transfer-matrix formahsm an exact description of the thermodynamic functions is obtained. The 
ground state phase diagrams for all regions of parameters values containing phases corresponding 
to the magnetization plateaus at M = 0, 1/4 and 1/2 have been obtained. Exact formulas for bond 
distortions for various ground states are presented. A mechanism of magnetization plateau stabi- 
lization with doubling of unit cell is reported. This mechanism is related with the inequivalence of 
right and left interaction bonds for certain lattice sites. 
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I. INTRODUCTION. 

The role of lattice distortions in the behavior of mag- 
netic systems continues to be in the focus of intensive 
theoretical and experimental investigations during last 
decade. The lift of the ground state degeneracy in frus- 
trated magnets and magnetization plateaus stabilization 
mechanisms have been investigated intensively in var- 
ious lattice spin systems with spin-lattice coupling^"— . 
The concept of adiabatic phonos^ yielding the effective 
spin Hamiltonian with additional biquadratic interaction 
is one of the main approximations using to gain inside 
into the properties of strongly correlated spin systems, 
interacting with lattice vibrations. Thus, considering 
only magnetic properties of the system with spin-lattice 
interaction, one deals with the Heisenberg model with 
additional biquadratic terms with coupling constant de- 
pending on the parameters of spin-lattice coupling and 
spring constant of the bond. Also, biquadratic terms 
can arise in the spin Hamiltonian from the quadrupole- 
quadrupolc interactions. The spin-lattice coupling in 
quasi-onc-dimensional frustrated spin chain, or zigzag 
ladder, has been shown to generate a series of magne- 
tization plateaus at various rational values of magneti- 
zation in addition to one at M = 1/3 existing in pure 
spin system^. It was also demonstrated than spin-lattice 
coupling in zigzag ladder gives rise to a novel type on 
magnetic excitations carrying fractional spin^. In two- 
dimensional both classical and quantum spin models with 
spin-lattice coupling the enhancement of magnetization 
plateau stabilization as well as appearance of new ordered 
phases due to biquadratic interaction have been also re- 
ported for the Ji — J2 — Js-model^ as well as for Shastry- 
Sutherland lattice^ and pyrochlore antiferromagneti. 

In this paper we consider an exactly solvable one- 
dimensional spin model with Ising and Heisenberg bonds 
and spin-lattice interaction between the spins connected 
with Heisenberg bond. The geometry of the model corre- 
sponds to the system known as sawtooth chain. Sawtooth 
chain, or delta-chain, is one of the prototype examples of 
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FIG. 1: The sawtooth chain. Heisenberg bonds are marked 
by thick lines, dashed lines represent Ising bonds. Group of 
sites pictured into the frame correspond to one block. 



highly frustrated lattices. It has a structure of corner- 
sharing triangles (Fig.(IT])). Unlike the other frustrated 
system, sawtooth chain has peculiar properties such us 
exactly known dimerized ground state and elementary ex- 
citations of quantum soliton type^rJ^. Significance of the 
Heisenberg model on sawtooth chain is not limited only 
to academical interest. Magnetic lattices of few class of 
materials, such as delafossite YCu02.3^'^ and olivines 
with structure ZnL2S4(L=Er,Tm,Yb)i^ have been found 
to be of sawtooth chain type. In addition to that, recently 
much attention has been paid to the problem of localized 
magnon states or dispcrsionless excitation bands, which 
have been found in various frustrated spin and electron 
systems, particularly, in the Heisenberg and Hubbard 
models on sawtooth chair^i^"— . Calculating the possible 
ground state degeneracy corresponding to the localized 
magnon (or electron in case of Hubbard model) states one 
can describe low-temperature thermodynamics near sat- 
uration field for the corresponding frustrated systemi^i^S. 
However, exact description of thermodynamic properties 
for sawtooth chain as well as for many other strongly 
correlated lattice model is still an open issue. It is worth 
mentioning two sophisticated methods which allow one to 
construct thermodynamic functions of integrable model 
- thermodynamic Bethe ansatz (TBA)^ and quantum 
transfer-matrix method (QTM)^, however, applicability 
of these methods is limited to a very narrow class of in- 
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FIG. 2; Local spin configurations corresponding to (from up 
to bottom) \AF2), \AF?,), \AF\), \SM) and \F\) ground 
states. Ovals denote \Slot ~ 0,+) singlet state of the pair 
of spins connected by XX Z bond. 



tcgrable systems, such as XXZ-Heisenberg chain with 
homogeneous couplings, Hubbard chain, e.t.c. Despite 
great successes in exact describing of thermodynamic 
functions for integrable models by TEA and especially 
by QTM2^i^, for many other physically and princi- 
pally important low-dimensional strongly correlated lat- 
tice models only laborious numerical calculations pro- 
vide more or less reliable results for finite T thermody- 
namics. Recently, many papers have been devoted to 
exact solution of the low-dimensional lattice spin mod- 
els with mixed Ising and Heisenberg bonds^i"— or just 
to pure Ising counterparts of known Heisenberg models 
with various one-dimensional topologies of bonds^"— . 
These exact solutions for modified models have much 
in common with the numerical and experimental re- 
sults obtained for their Heisenberg counterparts. For 
instance, an alternating spin-chain with ferromagnetic- 
ferromagnetic-antiferromagnetic(FFA) interactions has 
qualitatively the same magnetization curve with mag- 
netization plateau at M = 1/3 for Heisenberg^ and 
Ising spins^. For more complicated F-F-AF-AF chain 



it has been demonstrated that replacing all ferromag- 
netic bonds with Ising ones, holding the rest (antifer- 
romagnetic) bonds in original Heisenberg form does not 
lead to sufficient quantitative changes in the magnetiza- 
tion curve^i. The corresponding mixed Ising-Heisenberg 
F-F-AF-AF chain has been proposed in Ref. [23] as 
the model of magnetic structure of the compound Cu(3- 
Clpy)2(N3)2, where 3-Clpy indicates 3-Chloropyridine. 
These two examples are the simplest cases of various one- 
dimensional spin lattice models with clusters of spins in- 
teracting with Heisenberg interaction and interacting to 
each other only through Ising bonds, which have been in- 
vestigated intensively during last decade^!"— . Partition 
function and, thus, all thermodynamic functions for such 
models can be obtained analytically. All these facts allow 
one to consider these models as an approximation for the 
underlying quantum spin systems. In Rcf. [3^ the corre- 
sponding approximate model has been considered for the 
sawtooth chain. Sawtooth chain with Ising and Heisen- 
berg bonds considered there has four S" = 1/2 spins into 
the unit cell, thus, according to Oshikawa-Yamanaka- 
Afficck (OYA) criterion^ it can exhibit magnetization 
plateaus at M = 0,1/4 and 1/2 what has been found 
in Ref.jsl]. Plateau at M = 1/4 appears due to dif- 
ference into left and right interacting bonds for certain 
sites what, it its turn, results in unit cell doubling in the 
ground state corresponding to M = 1/4. In this paper 
we perform further investigations of this model comple- 
menting it with spin-lattice coupling and giving detailed 
analysis of all regions of the coupling constants values. 
The paper is organized as follows: in the Second sec- 
tion the Hamiltonian formulation of the model and its 
exact solution with classical transfer matrix methods is 
presented, the zero temperature phase diagrams are ob- 
tained in the third section, the fourth section contains 
an analysis of the average site displacements for differ- 
ent phases, concluding remarks are presented in the last 
section. 



II. HAMILTONIAN AND THERMODYNAMIC 
FUNCTIONS 



We start with considering the sawtooth chain with 
XX Z Heisenberg bond for the left pair of spins on every 
second triangle, while all other bonds are taken to be of 
Ising type. For the sake of briefness, hereafter, we will 
refer to these bonds as to quantum ones. Beside, each 
quantum bond is assumed to include elastic energy and 
displacement depended exchange constant (See Fig. [1]). 
Under these assumptions the Hamiltonian can be decom- 
posed into the sum of block Hamiltonians, commuting to 
each other: 



3 



^ / /I 



i=l 



Jn) SI, -{H- Ja,+i)S!^ + (1) 



r 



where N is the number of two-triangle blocks in the chain, 
Ji and J are coupling constants in the basement and be- 
tween the side sites and basement respectively, pi is the 
distortion of quantum bond on i-th block, A is dimen- 
sionless spin-lattice coupling, K is the spring constant 
and H is magnetic field assumed to point in z direction. 
^'i{i2) stand for a component of 5* = 1/2 first(second) 



spin operator of th block and arrangement of Ising a 
and T spins is depicted in Fig. ([T|). Due to commutating 
block structure of the Hamiltonian the partition function 
of the model can be represented in a way, suitable for ap- 
plying the transfer-matrix formalism. Namely, one can 
expand the exponential and obtain the product of the 
terms, with each one corresponding to one block: 



N 



^ exp (-[3 (Ui + JcTiTi ^ ^ yi + \ i^i + '^J+i) 



(2) 



r 



Here the sum is over all values of Ising variables and Sps 
stands for the trace over all states of quantum spins, for 
each quantum bond there is also integral over all possible 
values of distortion. After implementing integration, one 
arrives at the effective Hamiltonian with additional bi- 



quadratic term, which is widely used in the investigations 
of spin systems with adiabatic coupling to phononsi"— . 
Then, trace over quantum states can be taken indepen- 
dently for each block, yielding 



J 



N 



Z - ^ n 71=^^P[^^P ^)] exp (^-p (^Ja,n - H [r, + i (a, + ) ) ) , (3) 



where effective Hamiltonian reads 



= J (S,iS.2)a - b (S.iS.2)a -{H-Ji (cT,; + - Jn) Sf, - {H - Jct,+i)SI, 



(4) 



where the effective constant of biquadratic interaction is 
given by the relation b = ■ known, biquadratic 

term in case of spin s ~ \/2 yields the constant and bilin- 



ear term. So, one should also have in mind the following 
form of the effective Hamiltonian: 
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JXY = A{J 



JXY {SfiSf2 + Sf^S^2) + JzS-iS-2 - {H ~ Ji {ffi 



fjj+i) - Jn) S-i -{H ~ Jai+i)SI. 



Jz = J 



(5) 



Thus, one actually obtains the Hamiltonian of the XXZ- 
model with peculiar parametrization of the exchange 
constants and anisotropy. Any crucial changes in the 
properties of the system like appearance of novel phases 
which do not appear at 6 = docs not expected, while 
spin-lattice coupling can lead to some non-trivial im- 
pact on the ground states phase diagrams properties. 



Then, the partition function of that type can be calcu- 
lated exactly^"— within transfer-matrix formalism. To 
proceed, one needs to calculate the trace for each block 
quantum spins, and then to sum over two values of in 
each block to give Eq. (jS]) the standard form correspond- 
ing to the transfer-matrix technique^ 



T,=±l/2 



where 



n (a,, a,+i|T,) = e-^^"(^-^-+^l^') = 2e^T^ (e'f'i cosh (p\ {2H - J {r, + <j,+, - J, (a, + 



ri = l 



^^Pi{AH+j) ^^^^ hyiJ - n) - ,/i {a, + a,+,)f + (^J + j j , (7) 



where A„ {ai,ai+i\Ti) are four eigenvalues of "Hj 



Ai,2 (CTi,CTi+i|Ti) = (4J - 6) ± ^ (J (tj + cr,+i) + Ji {a., + tTi+i)) T H, 
Id 2 



A3,4 a,+,\T,) = (6 (1 + 4A2) + 4J) T - ^0 - 'h + <^^+l)f + J + 
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I 

The corresponding eigenstates are two polarized vectors states with S^^^ = for A3 4: 
with total spin projection S^^f. = ±1 for Ai^2 and two 



(8) 



5 



(in>-7±i;t»: 



7± 



J (CTi+l - Ti) - Jl {a.i + (Ji+l) ± Y ( J (CTi+l - T,) - Jl (cr^ + fTi+l)) + A2 (j + 



A J 



(9) 
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Thus, for calculating the partitfon functfon of the model, 
Eq. ([ni , one should find eigenvalues of the transfer matrix 
T, 



^0 



^0 



(10) 



where entries are connected with the value of block par- 
tition function as follows Z+ = Z(l/2,l/2),Zo = 



J 



Z(l/2,-l/2),Zo = Z(-1/2,1/2),Z„ 
Z(— 1/2,— 1/2). Having all that one can easily 
write the expression for the free energy per block in the 
thermodynamic limit, when only the maximal eigenvalue 
survives 
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One can see that the first additive term accounts for the 
vibration contribution to free energy , while effects of in- 
teractions between spin degrees of freedom and lattice 
displacements are incorporated into the 6-dependance of 
the rest part of free energy. After that, all thermody- 
namic functions can be found exactly by taking corre- 
sponding derivatives of the free energy. For instance, for 
magnetization one obtains: 



M 



2 \dHy J, 



(12) 



lational symmetry, and and 1 /4 for the state with dou- 
bling of lattice period. Accordingly, the magnetization 
curves with magnetization plateaus at all aforementioned 
values of magnetization are expected. Let us start with 
antiferromagnetic spin configurations with M = 0. There 
are three antiferromagnetic spin configurations with spa- 
cial period equal to lattice period in the system. In the 
first one all quantum spins are polarized, while both of 
their adjacent a and r spins point in the opposite direc- 
tion: 



III. GROUND STATES AND T = PHASE 
DIAGRAMS 

Analyzing possible magnetic ground states of the 
model under consideration at T = one finds large vari- 
ety of spin configurations with a spacial period equal to 
the period of the chain. However, there is also another 
two ground states where translational symmetry is bro- 
ken up to the period of two blocks. All states can be also 
classified by the values of magnetization, possible values 
of which are 0, 1/2 and 1 for states with unbroken trans- 



N 



\AFl)^Y[\n)^^W^=i,n=i), 



(13) 



here e is the corresponding energy per one block. Another 
two antiferromagnetic ground states of the system are 
degenerated by energy. They correspond to \S^g^ = 0, 
eigenstate (Eq. (|9]) ) for each XAZ-bond and opposite 
orientation for two Ising spins in each block: 
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\AF2) = n \SLt = 0, +). (g) I a, =t, n =;), 
1=1 

N 

\AF3) = n = 0' l"^^ =t)' 



1 , 1 . . . o^ 1 . .2 . ^ / - 1,^ ^ 



eAFi2,3)^-^J-Y^b{l + 4A')--d{J-Jir + A^lJ+-b] , (14) 



I 

where the value of the coefScient 7 in the ground state quantum spins are in {Sf^^ = 0,+) state while for each 

I'^tot — 07+)j is different according to Eq. (|9]). There is block (Ti=Ti but their joint direction alternates from block 

also another antiferromagnetic ground state with broken to block: 
translational symmetry. In this ground state pairs of 



JV/2 

\AF^) ^ n = 0, +)2.-l (g) k2,-l =t, =t) (g) l^fot = 0, +)2. (g) k2, =1, =1) , 

i=l 



The system demonstrates also three ferrimagnetic 
phases with value of magnetization equal to 1/2: 



AT 



1^1) - n \^tot = 0, +), (g =t, r, =t), 



AT 



i=l 

1 1 

£F2 = -T^o~ -Ji - H, 
lb 2 

AT 



£f3 = -T^b+l-Ji - H. 
lb 2 



(16) 



Local spin configuration corresponding to several eration are presented in Fig. ([2]). Finally, let us mention 
non trivial ground sates which appear as the zero- the spin polarized, saturated state, where all spins are 
temperature ground states of the system under consid- 
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pointed along the external magnetic field 

N 

\SP)=Y[\n)^^\<y^ =t,r,=t), 

esP = J+^Ji-^b-2H. (17) 

Within some region of the values of system parameters, 
which is specified below, there is also another state with 
broken translational symmetry. Namely, this is an inter- 
mediate state between \AF2) or |Ai^3) and \F1) where 



the (T spins are ordered antiferromagnetically with re- 
spect to each other, while all r spins are pointed along 
the field. This state can be achieved from \AF2) or |Ai^3) 
by flipping every second basement a spin. We will refer to 
this state as a spin-modulated stare as here the spacial 
modulation of local spin polarization is occurred along 
the chain. Namely, local spin polarization of each left 
pair of spins on each triangle varies in the following way 

100010001000 The corresponding ground state and 

energy per one block are 



N/2 



\SM) = n \SLt = 0, +>2,-l (g) 1^2.-1 =t, r2.-l =t> (g) \S^ot = 0, +>2. (g) W2^ =i, Ta. =t>, 



esM = ~^b{l + iA') - ij- i|A ( J+ ^fe ) I - + A'^ ( J+ ^6 



1 

2^- 



(18) 



Another one ground state can be appeared in the saw- 
tooth chain with Ising and Heisenberg bonds without lat- 
tice distortions^. When b = and J = Ji the ground 
state of the system is a macroscopically twofold degener- 
ated frustrated state 

N 

1=1 

where jf;) stands for cither \ai ^'t^n =1) or \ai =i,Ti =t 
). Thus, each pair of cr — r spin in each block is frustrated 
and can freely pass from one possible state to another. 

A. Ground state phase diagram for J > 0, Ji > 

Let us first consider the case of purely antiferromag- 
netic couplings J > 0, Ji > 0. The value of the bi- 
quadratic interaction constant is assumed to be always 
positive, as 6 = '^^^ ,K > 0. At such values of pa- 
rameters the system exhibits four different ground states, 
\AF2) or |Ai^3) which are degenerated by energies, \F1), 
\SM), which can appear only at r] = Ji/J < 1 and \SP). 
When the magnetic field in turned off the ground state 
of the system at T = is always an antiferromagnetic 
one, corresponding to either \AF2) or |^F3), which are 
closed analogs of dimerized ground states of conventional 
quantum sawtooth chain^"— . However, in the exactly 
solvable model considered here the structure of ground 
state is more complicated. In contrast to the conven- 
tional case, instead of simple dimers, ^ (| ti) — I it))i 
here one finds {Sf^^. = 0,+) states alternating with pairs 



of oppositely directed Ising spins Eq. (fT4| . Then, with the 
magnitude of the magnetic field increasing, a sequence 
of quantum phase transitions takes place, which increase 
the magnetization of the system in stepwise manner. The 
corresponding phase diagrams in the (r^, /i)-plane are pre- 
sented in Fig. ([3]), the upper panel demonstrates phase 
diagram for the sawtooth chain with Ising and Heisenberg 
bond without lattice distortions, which was previously 
obtained in Ref . [ [IHl , the lower panel demonstrates ef- 
fects of biquadratic term originated by lattice distortions. 
Actually, this term enhances the magnetization plateau 
at M = 1/2, which corresponds to the ferrimagnetic 
phase \F1). The region of phase diagram corresponding 
to the spin-modulated phase almost is not affected by the 
value of 6* = fe/ J. This can be explained by the fact, that 
spin-modulated phase which corresponds to the magne- 
tization plateau at the value M = 1/4, appears due to 
very special structure of the eigenstate of quantum spins 
with Sfgf. = which is strongly affected by the values of 
adjacent Ising spins. Absence of the permutational sym- 
metry in the {S^^^ = 0,-1-) state and specific dependence 
of the corresponding eigenvalue on the values of adjacent 
Ising spins makes the state with alternating orientation 
of basement Ising spins stable for certain value of param- 
eters, namely for 77 < 1. Thus, this is an example of 
the novel mechanism of magnetization plateau formation 
which consists in the doubling of the spacial period of 
the ground state due to the difference in effective fields, 
hi and /i2 which act on z-components of two quantum 
spins in one block, namely hi = H — Jti — Ji{(Ji + (J^+l)^ 
h2 = H — JiUi+i. The equations of phase boundaries 
between four aforementioned phases are: 
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FIG. 3: Ground state T = phase diagrams in {rj — Ji/ J, h — 
H/J)-plane for antiferromagnetic coupling (J > 0, Ji > 0) 
for the sawtooth chain with Ising and Heisenberg bond with- 
out biquadratic term (upper panel); and for b/J = 6 (lower 
panel). The corresponding eigenstates for the lattice are 
shown in Fig. Q 



between \AFi2, 3)) and \SM) : h ^ ^ (^l - \D+ {9) \ - ^1 + 0% {9)^ + ^ {I - i^f + 0% {9), (20) 



between \AF{2,3)) and |F1) 



h = 



^1 + ^/a vf + Dl (9) + D\ [9)^ , 



between \SM) and |F1) : 



h^-[l + \D+ {9)\ + -^l + D\ {9)]-^ii^ + Dl (9), 



between |F1) and \SP) : h = I + ^ U + ^tj^ + Dl (9) j + -A^6 



where 



Di (0) = ( 1 + -( 



(21) 
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-T/J=0.001 
T/J=0.01 
T/J=0.05 



Jl/J=0.5 
b/J=3 



1 

H/J 



-T/J=0.001 
T/J=0.05 
T/J=0.1 



J■^lJ=^ 
b/J=6 



0.0 



H/J 



FIG. 4: The plots of magnetization processes for the anti- 
ferromagnetic region of interactions, (J > 0, Ji > 0), for 
rj — 0.5 (upper panel) and = 1 lower panel for several tem- 
peratures. One can see effect of thermal fluctuation on the 
spin- modulated phase, the plateau at M — 1/4 shrink down 
very rapidly with increasing the temperature, while plateau 
at M = 1/2 is much more stable. 



Two typical plots of magnetization processes correspond- 
ing to all antiferromagnetic for different temperatures 
coupling are presented in Fig. (|4]), upper panel shows the 
77 < 1 with three magnetization plateaus at M = 0, 1/4 
and 1/2, the lower panel corresponds to the rj > 1 case 
with only two plateaus at M = and M = 1/2. As one 
can see spin-modulated phase in much more sensitive to 
thermal fluctuations than the rest plateaus, because it is 
connected with the configuration of Ising spins which is 
less stable that the entangled state of quantum spins. 

Looking at Fig. ^ one can also notice the effect of bi- 
quadratic interaction on the length of plateaus. It almost 
does not affect the plateau at M = 1/4 as spin-modulated 
phase appears due to spin-flip of Ising spins in the base- 
ment, however, minor deformations of the region corre- 
sponding to the spin-modulated phase can be observed 



from the phase diagrams. In its turn, biquadratic inter- 
action enhances the plateau at M = 1/2 corresponding to 
|F1) phase as well as the plateau at M = corresponding 
to \AF{2, 3)) phase. The effects of biquadratic term or 
lattice distortions arc summarized in Fig. , where the 
phase diagram in the {9, /i)-plane is presented. One can 
see almost parallel arrangement of the four phases men- 
tioned above and one peculiar point at = — 2 where 
spin-modulated phase disappears. Negative values of 
b are irrelevant in the model of lattice distortions, as 
b ~ and A' > 0. However, formally one can con- 

sider negative values as well. This effect is obviously 
connected with the peculiar properties of Hamiltonian 
with biquadratic interaction H ~ Ji (S1S2) + J2 (8182)^ 
at J2 = 2Ji. At this values all four eigenvectors be- 
came separable, corresponding just to the standard basis 
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FIG. 5; The ground states phase diagram in [Q = b/J,h = 
H/ J)-plane demonstrating effect of biquadratic terms for 77 = 
0.3 and A = 0.3. 



I tt) 1 1 i\) 1 1 ti) ) I it) - which makes impossible reahza- 



tions of spin-modulated phase. 



B. Ground state phase diagram for J < 0, Ji > 

Let us consider now another region of the coupling 
constants values corresponding to ferromagnetic J. In 
contrast to the previous case here three different ground 
states are possible at = depending on the value 
of 7] = jjj-. It is interesting to note, that the only 

one phase with zero magnetization (antiferromagnetic) 
which appears in this case is |Af 4) configuration with 
broken translational symmetry. Depending of the values 
of = yjj- and A the ground state phase diagram can 
have three qualitatively different forms (Fig |6|) , mainly 
determined by positions of two special points in 77 axis. 
One of them, rji, determines the quantum phase transi- 
tion point between spin polarized and antiferromagnetic 
|Ai^4) phases in the absence of magnetic field. Another 
one, ri2 corresponds to the position of quantum triple 
point in the (77, ft,)-plane: 



3 I 2 



Dt (6*) = A-' ( 1 - -t 



2 



e + Ji + D'i (9) 



A46'2 



(22) 



There is also another special point in the 77-axis cor- 
responding to the quantum phase transition point be- 
tween antiferromagnetic and ferrimagnetic ground states 
at H ~ 0. This value is parameter independent and is 
equal to one. When 7; > 1 the ground state of the system 
is |F1). One can see that at positive values of 6 always 
771 < 772 < 1. Thus, when 771 > one obtains the ground 
state phase diagram with spin polarized phase for 77 < 771 , 
antiferromagnetic phase |A_F4) turning immediately to 
spin polarized phase under the action of magnetic field 
for 771 < 77 < 772, antiferromagnetic phase |y4.i^4) turning 
to the ferrimagnetic phase \F1) prior to spin polarized 



phase for 772 < 77 < 1 and ferrimagnetic phase turning 
immediately to spin polarized phase for 77 > 1. One can 
see the corresponding phase diagram for particular value 
6* = 1 in the upper panel of Fig. Once 771 became 

negative, then one will get another topology of ground 
state phase diagram presented in middle panel of Fig. 
(|6|) for 6 = 2.5. And, finally, when 772 becomes negative 
only two ground states for H = are possible, |Ai^4) 
and \F1) and there is no quantum triple point in the 
phase diagram (see lower panel of Fig. where = 3). 
The equations of phase boundaries between three phases 
occurring in the negative J case are 
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between \AFA) and |F1) 



between {AFA) and \SP) 



1 / A^d 



between |F1) and \SP) : h ^ -1 



1 



A2 



+ Jl + D^_ (0) , 



C. Ground state phase diagram for J > 0, Ji < and \SP) one obtains 



(23) 



(24) 



At the ferromagnetic value of the coupling constant 
between the basement spins Ji < 0, the system exhibits 
rather simple ground states phase diagram with almost 
parallel arrangement of three phases, |AF(2,3)), |F1) 
and \SP). The H = ground state is always antifer- 
romagnctic which with increase of the external magnetic 
filed magnitude turned to the ferrimagnetic phase |-F1) 
with further transition to spin polarized saturated phase. 
The corresponding phase diagram is presented on Fig. 
d?]). Accordingly, the phase boundaries are given by the 
same equation as of the purely antiferromagnetic case 
J > 0, Ji > with replacement rj by —rj. (See Eg. (1^01) V 
Varying parameter of biquadratic interaction b docs not 
lead to a crucial change in the general picture. 



> JA 

(Po) = 



(27) 



while for other ground states average distance between 
quantum spins is given by more complicated expres- 
sions depending of the orientation of the other spins of 
the block. For the ground states with unbroken block- 
translational symmetry one obtains 



po = (flpol*) 

n 



(28) 



JA 



2A 



V 



IV. AVERAGE DISPLACEMENT 

One can easily determine the expression for the equilib- 
rium point of the quantum bond distortion immediately 
minimizing the block Hamiltonian ([1]): 

J A 

Po = -^(SiS2)^. (25) 

Thus, in the equilibrium, the distance between sites with 
quantum spins is determined by their magnetic state. 
In order to calculate the bond displacement for various 
ground states one just needs to take quantum-mechanical 
average of this operator. One can also exploit thermody- 
namical identities to determine standard deviation of p, 
more precisely 

On the other hand, one can obtain zero temperature val- 
ues of distance between the sites with quantum spins just 
by calculating quantum mechanical averages of the oper- 
ator pq — jj ^f=i Pio for corresponding ground states of 
the system. For all ground states where pairs of quan- 
tum spins are in | tt) state, i.e. for |^F1), |F2), \F3) 



where the following notations are adopted 

''^A(/+i6)' -"^^ A(J+i6)' ^''^ 

and l^") stands for the one of the \ AF{2, 3)), |F1) ground 
states, and 7* is the coefficient from the IS^^^ = 0,-1-) 
(Eq. ©) state calculated for the corresponding values of 
<Jl(R) E^nd r, where crL{R) is the value of the lcft(right) 
sigma-spin surrounding quantum spin pairs in the given 
ground state. For the ground states with block doubling, 
|^F4) and \SP) the corresponding expression has the 
form 

(*l/5o|*) = ^(pi+P2), (30) 

where pi^2 £^rc calculated for the first(sccond) block in 
the two-block unit cell of the ground states configuration. 
Each Pi is given by Eq. ([^5]) with corresponding values of 
'^L(R) T. In order to compare zero-temperature re- 
sults with thermodynamically obtained expression form 
Eq. (|26p one needs also to calculate quantum-mechanical 
averages for the square of operator po . Then, one will ob- 
tain a quantity thermal average of which coincides with 

^ from Eq. (gB]). Defining ^0 = \J (^Ip^ol*) o^e will 
obtain for ground states without translational symmetry 
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breaking 



4K 



\ 



1 + 4A 



V 



and for |AF4) and \SP) phases 



(*|p\l*> 



where ^1(2) are calculated according to Eq. ([31]) for left 
and right triangles in the block. Below the list of dis- 
(Sl^ances between the site with quantum spins for different 
Tphases for zero temperature is presented 



1 + ( J(r - an) - Ji{(Jl + cfr) 



(32) 



|^F(2,3)): po 



\AFA): p, = - — 



\Fl) : po = -7T7 




\SM) : pi 
JA 



JA 

'4K 



2A \ , JA L ^ 1 

, ^ , p2 = -—{l + 2A), p, = - [l + A[l + ^= 



4a: 



1 + 4A 



(33) 



1 + 2A 



2A 



Thus, from Eqs. (|27| and (p3| one can conclude that 
for J > the interplay between elastic properties of 
the quantum bonds and magnetic behavior od spins con- 
nected by it yields the stretching of the bond length for 
those ground states in which S^^^ = 1 for each quantum 
bonds (I^Fl), |F2), \F3) and \SP)), while for the ground 
states with S^^t = (|Ai^(2,3)), \AF4}, \F1) and \SM)) 
antiferromagnetic coupling of spins connected by XXZ 
bond leads to shortening of bond length. Though, in 
later case the expressions for the average displacement 
are moro complicated and contain all interaction param- 
eters of the system, only the sign of J and A defines 
whether the equilibrium displacement will be positive or 
negative. Magnetic field dependence of the standard de- 
viation of the distance between lattice sites with quantum 
spins given by Eq. is presented in Fig. ^ . One can 
see low temperature step-like changes of ^ according to 
the connection between site distance and magnetic state 
of spin situated at them (Eq. ([25])). 



V. CONCLUDING REMARKS 

In this paper we presented a complete analysis of the 
magnetic properties of the sawtooth chain with Ising and 
Heisenberg bond and spin-lattice coupling for spins in- 
teracting with Heisenberg interaction. After integration 
over site displacements one deals with the additional bi- 
quadratic spin interaction. Due to special arrangement 
of Ising and Heisenberg bonds exact calculation of par- 
tition function of the system has been performed. The 
system exhibits large variety of ordered phases, among 
which are those with the doubling of the unit cell. Mag- 
netic behavior of the system is also rather rich. Depend- 
ing of the values of parameters, coupling constant in the 
basement Ji, coupling constant for the interaction be- 
tween spins on the top and basement spins J and effective 
coupling constant of additional biquadratic interaction 6, 
various magnetization curves are possible with magneti- 
zation plateaus at M = 0, 1/4 and 1/2. Let us mention, 
that exact diagonalization calculations for ordinary anti- 
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ferromagnetic sawtooth chain at J = 2 Ji reported in Ref. 
[2l| . revealed only one plateau at M ~ 1/2, which cor- 
responds to the so-called magnon crystal, an eigenstate 
with whole filling of all possible localized magnon states. 
In our case, the microscopic physical origin of this plateau 
state is completely different due to special structure of 
the interactions. However, at corresponding values of 
coupling constants (the left half of the phase diagram 
presented in Fig. (jSj in upper panel) the system with 
mixed Ising and Heisenberg bonds considered here ex- 
hibits magnetization curve with all three magnetization 
plateaus mentioned above. These features have obvious 
origin. The ordinary sawtooth chaiit^i has unit cell with 
two S" = 1/2 spins, thus, according to OYA criterion^, it 
can display only plateaus at M = and M = 1/2, while 
the system considered in this paper has four S* = 1/2 spin 
in the unit cell which can result in additional magneti- 
zation plateau at M = 1/4. However, OYA criterion^ 
only specify possible values of magnetization plateaus for 
given type on translational symmetry and given values of 
spin. For stabilization of the possible plateaus additional 
physical mechanisms are required. In the system con- 
sidered here a special attention should be paid to the 
magnetization plateau at M = 1/4. The origin of the 
corresponding eigenstate \SM) is in the nonequivalence 
of the left and right a spin for each pair of quantum spins, 
provided later are in the jS'^^o^ =0,-1-) state in which Si 
and S2 in their turn are nonequivalent. Thus, at a certain 
value of external magnetic field magnitude, the eigenstate 
with doubling of the unit cell becomes stable, in which 
cr spins in the left and right blocks of the doubled unit 
cell are in different states. To our knowledge, this is the 
novel mechanism of magnetization plateau stabilization 
inherent in the systems with mixed Ising and Heisenberg 
bonds. The same arguments are valid for the antiferro- 
magnetic state with doubling unit cell \AFA). Generally 
speaking, the same mechanism can be presented in the 
purely quantum models with inhomogeneous interaction. 
For sawtooth chain, one can consider the model where 



coupling constants in the basement, on the left leg and 
on the right leg of the triangles, are all different, or con- 
sider the model with alternating interaction in the base- 
ment. Probably, for some values of coupling constant 
the magnetic behavior of such a sawtooth chain will have 
much in common with that shown here for model with 
Ising and Heisenberg bonds, particularly, magnetization 
plateau at M = 1/4. This issue requires further investi- 
gation in order to clarify the deep connections between 
magnetic and thermodynamic properties of quantum spin 
models and their Ising-Heisenberg counterparts. There is 
surprisingly good correspondence between magnetization 
curves of quantum F-F-AF-AF alternating chain and the 
same chain with ferromagnetic bonds changed with the 
Ising ones obtained in Ref. [13]. This, as well as other 
rcsultsSZr— , allows one to consider the change of inter- 
action bonds with Ising ones for providing exact solv- 
ability, as an approximate methods in theory of strongly 
correlated spin lattice models. An especially good agree- 
ment can be achieved when only ferromagnetic bonds are 
changed with Ising ones, because the ground state of two 
fcrromagnctically interacting spins is the same for Ising 
and for Heisenberg interaction. While, for antiferromag- 
netically interacting spins quantum ground state is spin 
singlet which is an entangled state which has no direct 
analogies in case of Ising interaction. Thus, changing an- 
tiferromagnetic bond with Ising one leads to crucial loose 
of important physical properties, while for ferromagnetic 
coupling the physical difference is not so pronounced. 
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FIG. 6: The ground states phase diagram for J — —1, A = 1.2 
and 9 = l(upper panel), 9 — 2.5 (middle panel) and 9 — 3 
(lower panel). Here h — H/\J\,9 = h/\J\. 
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FIG. 7: The ground states phase diagram for Ji — —1, J = 1, 
e = 1 and A = 1. Here v = \ Ji \/J, h = H/J. 
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FIG. 8: Standard deviation of the distance between sites with 
quantum spins vs. external magnetic field magnitude for three 
different temperatures and J — 1, Ji = 0.5, 6 = 0.5 and 
A = 0.3. Regions of different ground states corresponding 
to T = phase diagram are separated from each other by 
horizontal thin lines. 



